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ABSTRACT 

Extra-solar planets direct imaging is now a reality with the deployment and commis¬ 
sioning of the first generation of specialized ground-based instruments (GPI, SPHERE, 
P1640 and SCExAO). These systems allow of planets 10 7 times fainter than their host 
star. For space-based missions (EXCEDE, EXO-C, EXO-S, WFIRST), various teams 
have demonstrated laboratory contrasts reaching 10“ 10 within a few diffraction limits 
from the star. However, all of these current and future systems are designed to detect 
faint planets around a single host star or unresolved multiples, while most non M-dwarf 
stars such as Alpha Centauri belong to multi-star systems. Direct imaging around bi- 
naries/multiple systems at a level of contrast allowing Earth-like planet detection is 
challenging because the region of interest is contaminated by the hosts star companion 
as well as the host Generally, the light leakage is caused by both diffraction and aberra¬ 
tions in the system. Moreover, the region of interest usually falls outside the correcting 
zone of the deformable mirror (DM) for the companion. Until now, it has been thought 
that removing the light of a companion star is too challenging, leading to the exclusion 
of binary systems from target lists of direct imaging coronographic missions. 

In this paper, we will show different new techniques for high-contrast imaging of plan¬ 
ets around multi-star systems and detail the Super-Nyquist Wavefront Control (SNWC) 
method, which allows to control wavefront errors beyond nominal control region of the 
DM. Using the SNWC we reached contrasts around 5 x 1CU 9 in a 10% bandwidth. 

Subject headings: exoplanets, Alpha Cen, double stars, extended disks, wavefront con¬ 
trol, MEMS 


1. Introduction 


The exoplanets field is rapidly expanding with the success of the Kepler mission ((Burke 
et al. 2014) and references therein) and the emergence of direct imaging ground based instruments 
(GPI ( Macintosh et ah]|2014 ), SPHERE ( Beuzit et al.|[2008 ), SCExAO ( Guyon et al.| 2010), P1640 
(Hinkl ey et al.| 2008). One of the most exciting prospects of future telescopes is finding other Earths 
analogues in our galaxy or solar neighborhood and more ambitiously detect life on them. The Kepler 
space telescope has already revealed that roughly 22% of stars have planets between 1 and 2 Earth 
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radii in their habitable zone (Batalha 2014). However, this mission does not perform spectral 
characterization of these targets. Direct imaging combined with spectroscopic characterization of 
exo-Earths would allow us to determine the chemical composition of the planet’s atmosphere, and 
constrain the presence of oxygen, water and other elements necessary for life. Over the past decade, 
there have been more than a dozen direct imaging mission studies for space-bases telescopes. Figure 
[l] shows a few representative missions into which a lot of these have evolved. 



aCen imager 

Belikov et al 

(25cm) 


EXCEDE 

Schneider et al. 

(0.7m) 


Exo-C 

Stapelfeldt etal. 

(1.5m) 


Exo-S 

Seager et al. 

(1.1m w/ starshade) 


AFTA 

(2.4m) 


ATLAST 

( 8 - 16 m) 


Fig. 1.— Mission concepts compatible with the methods proposed here. For these missions, no 
hardware changes are required as long as they have a deformable mirror. 


An Earth-like planet orbiting the habitable zone of a sun-like star would be 10 billion times 
dimmer than the star. Diffraction created by the telescope aperture as well as aberrations is 
several orders of magnitude brighter than the planet, making its detection very challenging. Several 


starlight suppression systems (Guyon et al. 


2010 


Kern et al. 2013) have been demonstrated 10 


-9 


raw contrast or better in the laboratory. 


These systems employ high performance coronagraphs to suppress the star diffraction created 
by the telescope aperture and efficient wavefront control systems based on a deformable mirror 
(DM) to remove residual starlight leak (speckles) created by the imperfections of telescope optics. 
However, their designs have thus far been mostly limited to single-star systems and to planets or 
disks found within the DM control zone, where speckles can be corrected by conventional means. 
Binary star systems are good candidates to search for planets since they are more common than 
single stars. Good examples of such binary systems are Alpha Centauri (a Cen) and Sirius. Cur¬ 
rently, these multi-star systems are typically excluded from mission target lists because there is 
no technical approach that can manage the technical challenges associated with double-star (or 
multi-star) high-contrast imaging. The three main challenges are: 

- The first challenge is that often the multi-star separation is typically beyond the half-Nyquist 
frequency of the deformable mirror. In this paper, we propose a new method called Super-Nyquist 
Wavefront Control (SNWC) that uses a mild grating (or an existing pattern, print-through, com¬ 
monly found on many Deformable Mirrors left over from their manufacturing process) to effectively 
alias low-spatial frequency modes of the DM into higher frequencies, enabling the DM to remove 
speckles well beyond the DM’s Nyquist frequency. In effect, aliasing is used as a feature rather 
than a bug. 

- The second challenge is to separate and independently remove overlapping speckles from multiple 
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stars. The solution is called Multi-Star Wavefront Control (MSWC). We solve this challenge by 
selecting a region of interest such that different DM modes are used for the two different stars (at 
the expense of reducing the discovery region). 

- The third challenge is to create a dark zone for a multi-star system where the planet angular 
separation star requires SNWC to control diffraction in the discovery region, and the companion 
star creates diffraction and aberration light in the same discovery zone. In this case we combined 
the SNWC and MSWC together, and call this a technique Multi-Star Super Nyquist Wavefront 
Control (MSSNWC). 

Each of these techniques serve specific science case that we will present in section [2} In this 
paper, we focus on the Super Nyquist algorithm. After discussing the challenges of observing multi¬ 
star systems in section |3j we present the theoretical background of the method in section |4j We also 
briefly explain the multi-star wavefront correction. However the details of the MSWC is beyond 
the scope of the paper. Finally, we show simulation results in section [5] both in monochromatic and 
polychromatic light. 


2. Science motivation 

2.1. Searching for planets around multi-star systems 

The science cases addressed by several exoplanet detection missions such as the ones shown 
in Figure [l] may result in great leaps in our understanding of warm disks, exoplanet diversity, 
dynamics, and atmospheres. They also will deliver a census of exoplanets around nearby stars, 
and (for some missions) the detection and spectral characterization of Earth- like planets in the 
habitable zone. However, none of these missions are planning to image multi-star systems with 
current technology, except systems for which the leak and glare of the companion star is negligible. 

The method presented in this paper will greatly multiply the science yield of these missions 
since it enables direct imaging of planetary systems around multi-star systems, without additional 
hardware changes or costs. Enabling the study of multi-star systems is very important because the 
majority of K-type and earlier stars are in multi-star systems. In particular, 5 out of the 7 star 
systems within 4 parsecs containing K- or earlier type stars are multiples (a Cen, Sirius, Procyon, 
61 Cyg, e Ind), and only two are singles (e Eri, r Cet). While it is true that there are many more 
nearby M-dwarfs and most of them are isolated, the direct imaging of M dwarfs is arguably better 
done from the ground with Extremely Large Telescopes (ELTs) because: (a) they are dimmer and 
require larger apertures; (b) their planetary systems are closer to the star and require the angular 
resolution of larger apertures; Conversely, the study of K- and earlier type stars arguably favors 
space-based missions, for which the smaller apertures are less concerning, but require much deeper 
contrasts only possible from space. Therefore, most of the stars best suited for space missions are in 
fact in multi-star systems and SNWC combined or not with MSWC promise to greatly increase the 
science yield of almost any direct imaging mission, as well as enable the study of a whole new class 
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of star systems, namely multiple stars. One case deserves special attention. SNWC and MSWC 
enable high contrast imaging around our nearest-neighbor star, aCen. High-contrast imaging with 
an inner working angles of 2 X/d are now realistic, a 30cm telescope should be sufficient to image 
the habitable zone at high-contrast. The aCen system represents a particularly favorable target for 
planet detection missions, but it has been excluded from current missions’ target lists because of 
its multiplicity. Stars of comparable proximity to a Cen are all very dim, and stars of comparable 
brightness are about three times farther away. In particular, the next closest star earlier than M- 
type (e Eri) is 2.4 times as far, and is known to have a thick disk that may interfere with detection 
of small planets. The next star of comparable proximity to a Cen is Barnards star, which is 1.4 
times farther, has a much dimmer magnitude (M=10), and has a habitable zone only 30 mas wide, 
requiring at least a 4m aperture to even resolve it (Figure [2]). 



Fig. 2.— Simulation showing an Earth twin at maximum separation around every real nearby star. 
On this contrast vs separation angle graph, the circle size and color represent star size and type. 
aCen is 3 times easier than any other star by almost any metric, except for the fact that it is a 
binary. The vertical line also shows the particular case of where 2 X/d would be placed for a 1.5m 
aperture telescope at 550nm. 


Another reason why a Cen is an attractive target is because recent estimates of r/@, from Kepler 
have been on the order of 10% per octave of semi-major axis and per octave of planet size, leading 
to integrated numbers of 10- 22% under typical assumptions of habitable zone size and habitable 
planet size range (Batalha 2014). Therefore a Cen has about a 20-40% chance of harboring an 
exo-Earth around either the A or B star. A mission using the techniques proposed here, may be 


the first to detect and spectrally characterize an Earth twin, if one exists around a Cen (Dumusque 




















2012 ). 


et al. 


2.2. Exoplanetary debris systems 


SNWC allows another specific science case. When using SNWC, the outer working angle 
(OWA) and thus the field of view (FoV) of most missions in Figure [l] is limited by the DM half 
Nyquist frequency. As an example, for 2.4m telescope with a 48x48 DM such as AFTA config¬ 
uration, this implies an outer working angle of about 1” (depending on wavelength). Thus, the 
disk around e Eri can only be seen in high contrast out to ~4AU. For the case of aCen, AFTA 
would not be able to observe beyond 1AU even if the second star was not there. The SNWC 
method resolves this outer working angle limitation since it enables the extension of the dark zone 
for single-star systems past the DM Nyquist limit. SNWC does not increase the size of the dark 
zone, but it does allow a user to shift it to arbitrary locations. By stitching separately acquired 
sub- and super-Nyquist dark regions, imaging of arbitrarily large disks is enabled. For instance, 
simulations Q show that HR4796A, the disk is truncated with the current capabilities of the AFTA 
coronagraph instrument (Figure [3]). 




Fig. 3.— Simulations (courtesy of Thomas P. Greene, private communication) 
for the example of HR4796A, the disk is truncated with the current capabilities 
instrument.The simulations were done using the Zodipic package, a general-purpose 
for optically-thin disks (Kuchner 2012). 


showing that 
of the AFTA 
modeling tool 


3. Challenges of multi-star and extended disks observing 

As mentioned in the introduction, coronagraphy is challenging in the context of multiple-star 
systems because even in the case of external coronagraphy, the companion star leaks light into the 
region of interest around the target. Moreover, the separations between the two stars is such that 


1 courtesy of Thomas P. Greene, private communication 
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it is outside the correctable field of view (FoV) set by the number of degrees of freedom of the 
deformable mirror. In the context of extended disks, the limitation is the limited FoV. 
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Fig. 4.— Schematic showing the limitations of the multi-star coronagraphy. The image shows 
overlapping speckle fields from the two companions (A and B). The black squares show the standard 
control zone of the deformable mirror with respect to each star, which we refer to as the ”sub- 
Nyquist region” (of that star) and the rectangle called DZ is the dark zone region, where we want 
to detect a planet. The mean and median contrast of this dark zone are respectively 5.4 xlO -5 
and 3.49 x!0“ 5 
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Light leakage from the companion star(s) When observing a target belonging to a binary 
system, the amount of light leaking from the off-axis companion (hereafter referred to as the 
companion) reduces the Signal to Noise Ratio (SNR) of the candidate planet around the on-axis 
target (hereafter referred to as the target) to a level that the planet candidate might become 
undetectable. Two fundamentally different effects contribute to creating this light leak: diffraction 
and optical aberrations. These effects are complicated by the fact that the beams coming from the 
two stars are incoherent from each other. 



Distance from the star in I/d Distance from the star in I/d 


Fig. 5.— Light diffracted at 100 X/D in monochromatic light, without aberrations (left) and with 
25nm mis aberrations (right). The size of the region of interest is 4x8 X/D. The top views are 
larger field of views showing the companion star. The bottom is a zoom on the region of interest 
(close to the suppressed on axis star). The median contrast intensity of the order of 4x 1CF 8 
without aberrations and 4x 10 - ' with 25nm rms of aberrations. The contrast is shown in a perfect 
scenario (meaning totally removing the light from the on-axis target). 


To quantify the light leakage, we consider the case of searching for a planet or a disk located 


























at 100 A/d away from a target star. This could be a binary separated by 10” observed through a 
1.5m telescope at a wavelength of 770 nm or a single star observed with a large telescope at shorter 
wavelength.Figure [5] shows the light diffracted by the companion in the particular scenario where 
the on-axis star is totally suppressed (valid in the context of EXO-S for instance, for which the star 
shade practically removes all the light from the target star). The simulated contrast is the intensity 
at a given position in the image normalized by the peak intensity of the target. The simulation was 
done in the context of no aberrations (left image) and with 25nm rms of aberrations with a power 
law of /~ 2 (right image). The figures also show the dark zone region (DZ) as a box. We use these 
regions as a reference for the initial median contrast before corrections. 


mean: 1.17e-07, median: 1.18e-07,4.0 - 8.0 TO 
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Fig. 6.— Light diffracted past 100 X/D away from a star in polychromatic light (10% bandwidth), 
with lOnm rms aberrations. The size of the region of interest is 4x4 X/D. The median contrast 
intensity is of the order of le~ 7 with lOnm rms of aberrations, masking the light coming from a 
potential planet. 


The amount of light originating from the companion depends on the aberrations amplitude 
and their spatial frequency distribution in the system. In monochromatic light and in presence of 
no aberrations, the amount of light diffracted in the DZ is of the order of 4e -8 . With 25nm rms 
of aberration, the contrast worsens to 4e - '. Figure [6] shows the polychromatic case, for a 10% 
bandwidth around the central wavelength of 550nm. The median contrast intensity is of the order 
of 2e -8 without aberrations and 2e -7 with lOnm rms of aberrations. In monochromatic light the 
diffraction rings are blurred together and therefore not shown on Figure [6j 

Until now the standard approach to control diffracted light from the second star in a binary 






9 


system has been to design a coronagraph to block both stars. A coronagraph is only useful if the 
diffraction dominates over aberrations or otherwise cannot be removed by the wavefront control 
system. As a general rule, for typical mirror aberrations, diffraction dominates close to the star but 
aberrations dominate far from the star. All coronagraphs suppress diffraction, but they do not help 
with aberrations. A wavefront control (WC) system (or equivalent) is required to suppress optical 
aberrations, and in addition to suppressing aberrations, is often capable of suppressing diffraction 
by at least an order of magnitude. Therefore, a coronagraph is neither sufficient nor necessary to 
suppress the leak of the off-axis star, while a wavefront control system is both necessary and often 
sufficient to suppress both the aberrations and diffraction. 


Sub-Nyquist region of the DM The nominal region over which we can create a dark zone 
is defined by the number of actuators on the deformable mirror. Figure [7] illustrates this issue. 
The square region around the central star is the controllable region and is limited by the Nyquist 
frequency /at, which corresponds to an outer working sky angle of to N ac t x A/2 D. It has been 
previously believed that the deformable mirror can not control any speckles past this Nyquist- 
limited outer working angle. 



Fig. 7.— Simulation of a dark zone with a discovery region limited by the number of actuators on 
the deformable mirror. This is a case where only phase is corrected leading to a symmetrical dark 
zone. The coronagraph blocks most of the light before an image plan wavefront control algorithm 
further removes speckles. 
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4. Theory of Super-Nyquist Multi-Star Wavefront Control 
4.1. Super-Nyquist Wavefront Control 

We propose here to develop a solution based on wavefront control rather than using a new 
coronagraph design. Our solution to the super-Nyquiest limitation uses a mild grating (or an 
existing pattern commonly found on many DMs left over from their manufacturing process) to 
remove light coming from the off-axis star. The grating creates a controlled diffracted image of the 
off-axis star close to the on-axis star, placing the diffracted off-axis image within the DM control 
region, enabling it to remove the speckles from both stars (which can be done independently at the 
cost of halving the size of each control region). This solution uses the DM in a regime outside its 
nominal control range, where the DM experiences what is typically a harmful side effect: spatial 
aliasing. We use aliasing as a feature. 

In this section we present the theoretical foundation for Super-Nyquist Wavefront Control and 
explore its theoretical capabilities and limitations. We also discuss how to design DM parameters 
for a particular application. Before we begin, we introduce a simplified formalism for general 
coronagraphic wavefront control. 

Consider some arbitrary coronagraph with a deformable mirror and a science focal plane. Let 
Edm{x ) and Ef(£) be the electric fields in those two planes, where x and £ are 2-vectors representing 
the normalized 2D coordinates in units of pupil size D and sky angle A /D, respectively. Because 
the coronagraph is a passive linear system, the relationship between these two fields is given by 
some linear operator C: 

Ef = C{Ejjm} (1) 

A change in the deformable mirror setting creates a change A Edm in the DM field and a 
corresponding change A Ef = C{AEdm} in the focal plane. 

For purposes of our simplified treatment, we will assume that most coronagraphs can be 
approximated as follows: 

C{E dm (x)}(0 = F{A{x)E dm {x)}T{0 (2) 

where T is the Fourier Transform (with our normalization of x and £, the Fraunhofer integral reduces 
to the Fourier Transform to within a constant phase factor), A{x) represents the coronagraph 
aperture (or more generally, the cumulative effect of all apertures projected or propagated onto the 
DM), and T(£) represents the coronagraphic throughput as a function of sky angle, normalized to 
a maximum value of 1. Typically A{x) = 1 across a large portion of the DM, and T(£) = 1 almost 
everywhere in the focal plane except the small blind spot with radius of a few X/D corresponding to 
where the star is suppressed by the coronagraph. In what follows, we make the following simplifying 
assumptions: (a) we only consider the regions outside the coronagraphic blind spot, which allows 
us to set T(£) = 1, and (b) the aperture A(x) is binary valued (0 or 1), and is 1 across most of 
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the aperture. Note that with these assumptions and our normalizations, Eq.Q will be in units of 
contrast. With these assumptions and simplifications, we get: 

A E f = E{AE DM } (3) 

The key idea behind SNWC is that a DM can diffract some light in A Ef beyond the sub- 
Nyquist region into super-Nyquist regions, thus enabling control of those regions. Here, by ”sub- 
Nyquist” we mean the region within the spatial half-Nyquist frequency of the DM (corresponding 
to Oth diffraction order), and by ’’super-Nyquist” we mean regions beyond the spatial half-Nyquist 
frequency of the DM (corresponding to 1st and higher diffraction orders). For example, for a 32x32 
actuator DM, the sub-Nyquist region would be within 16A/D outer working angle. Conventional 
high contrast wavefront control is limited to this sub-Nyquist region only. However, because DMs 
diffract and modulate light in super-Nyquist regions, stellar leak in those regions can in principle 
be controlled in essentially the same way as in the sub-Nyquist region. The main difference is 
that the contrast level of DM modulation will generally be weaker in super-Nyquist regions than in 
sub-Nyquist because the amount of light diffracted into those regions by the DM is generally small. 
We will quantify and show how to mitigate this effect as part of the analysis in this section. 

Light can be diffracted either due to, (a) the periodic actuator nature of the DM, including 
any print-through pattern, (b) external gratings (which may not share the actuator periodicity). 
We treat those separately in the two subsections below. 


4.2. SNWC with diffraction caused by the DM influence function or print-through 

pattern 


The DM field is given by: 

E dm (x ) = A(x)e^ aM ^ (4) 

= A(x)[l + i<f>DM(x) + o((/) D m(x) 2 )} 

where 4>dm(x ) is the phase imparted to the electric field by the DM in radians. The first term of 
the above equation corresponds to the on-axis PSF, while the rest is the contribution is due to the 
DM which we will call A Edm- For small DM modulations, we can assume that only the leading 
term is significant and use the influence function model: 

A E dm (x) = i(j)DM(x) 

N 

= i ^2 a n f{x — nd) (5) 

71=1 

where / is the DM influence function, d is the spacing of actuators on the DM, and a n (for n = 1...N) 
are the DM actuator coefficients. We also dropped A(x) for simplicity and will simply assume the 
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equivalent condition that 4>dm(x ) = 0 outside the DM aperture. We will adopt the convention of 
/ being normalized to unity maximum value, which implies that a n are in units of radians. We can 
re-express this equation as a convolution: 


A E dm (x) = ia* f 


where a 


N 

52 a n 5{x — 

71=1 


nd). This convolution can be seen on the left column of Figure 


8 


. We now 


proceed to analyze the behavior in the focal plane by computing the electric field change created 
by the DM: 


A Ef = J-’jA Edm} 

= iF{a}F{f} (6) 

Computing the intensities (which with our normalized units will be in units of contrast) gives: 

\AE f \ 2 = \F{a}\ 2 \F{f}\ 2 (7) 

E{a] is a periodic function (see Figure [8] , top row), because the periodic sampling of a by delta 
functions makes a n the Fourier Series of Fa. Its different periods correspond to the sub-Nyquist 
region periodically copied, or aliased, into super-Nyquist regions. In other words, if influence 
functions were delta functions (a non-physical hypothetical scenario), any speckle modulations by 
the DM in the sub-Nyquist region would be perfectly periodically repeated (aliased) into all super- 
Nyquist regions. In the case of realistic influence functions, the super-Nyquist regions will get 
attenuated (as we show below), imposing a coupling between the width of the influence function, 
DM stroke and maximum correctable contrast of errors in super-Nyquist regions. This coupling is 
not simple to characterize in the general case, but becomes very simple if we characterize it in a 
statistical sense. We can treat the DM actuator coefficients a n as independent random variables 
with standard deviation \J(a^) radians (where without loss of generality we assume a 0 mean). 
This standard deviation is essentially a measure of DM stroke. For example, in a 32x32 actuator 
DM, the peak stroke will typically be 3x this amount. 

Random variations of the DM lead to random variations of the focal plane field Ej (” speckles”). 
The power spectral density of these variations is: 


(\AE f \ 2 ) = (\F{a}\ 2 )\F{f}\ 2 (8) 

N 

= Y,( a n)\F{5(x-nd)}\ 2 \F{f}\ 2 

n=l 

= N(a 2 n )\F{f}\ 2 


This power spectral density is essentially a measure of speckle contrast (as a function of position 
in the image plane) correctable by a DM with stroke yj (a^). As expected, higher strokes lead to 
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higher energy or contrast in the DM correction field and therefore the ability to correct for higher 
levels of error. Because the vast majority of the DMs are capable of strokes of several radians 
(yj (aj*) >1), and because our linearity assumption breaks down for {yj (a^) > 1 ), we can assume 
that for the vast majority of cases, y/ (a^) = 1 , simplifying the above to: 

(|A£/| 2 ) = N\T{f}\ 2 (9) 


which states that the Fourier Transform of the influence function times the number of DM actuators 
N directly gives a measure of the contrast correctable by the DM, in sub- as well as super-Nyquist 
regions. 


* 




ae dm (P hase > 
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Fig. 8.— DM print-through case. The left column shows the pupil plane and the right column 
the image plane. The top row shows the actuator coefficients a and its Fourier Transform (Fourier 
series of a n ). The second row shows realistic influence functions and the bottom row shows the 
convolution or multiplication of the top two rows. 


Figure [8] shows practical applications and implications of this equation. The left column shows 
the pupil plane and the right column the image plane. The top row shows the actuator coefficients a 
and its Fourier Transform (Fourier series of a n ). This would be the situation if the influence function 
was a delta function (unphysical). In this case, the correction from a DM in the focal plane is a 
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periodic function repeating in each Nyquist-sized region. Controlling or modulating speckles in 
the central sub-Nyquist region would be perfectly repeated in all super-Nyquist regions, enabling 
independent suppression of starlight in any super-Nyquist region in exactly the same fashion as in 
the sub-Nyquist region. The second row shows realistic influence functions and the bottom row 
shows the convolution or multiplication of the top two rows. Note that the middle right plot also 
represents the contrast of error correctable in super-Nyquist regions per Q. We now consider four 
representative cases of influence functions. 

Typical DM influence functions in continuous-sheet DMs (Figure [8j middle left, blue) are 
roughly the size of actuator spacing. The corresponding (|A£j| 2 ) (Figure [8j middle right, blue) 
tells us that there is only significant correction in the sub-Nyquist region while super-Nyquist regions 
are attenuated too much to be controllable at any reasonable contrast levels. Roughly speaking, 
because such an influence function is a factor of N smaller than the DM itself (in area for 2D), © 
tells us that the width (area) of (|A£/| 2 ) is NX/D ( N(X/D ) 2 for 2D) and the mean value is 1/N. 
This simply states the well-known fact that conventional DMs are capable of suppressing starlight 
leak in the sub-Nyquist region, and that the mean contrast of that leak cannot be higher than 
1/A, averaged across the sub-Nyquist region, by energy conservation. For example, in a 32x32 
DM (N = 1024), the sub-Nyquist region is 1024(A/D) 2 in area (32 x 32 A /D) and the controllable 
level of starlight leak is 10 -3 contrast averaged across the sub-Nyquist region. For such a DM, 
super-Nyquist control requires a grating or print-through pattern on the DM (covered in the next 
subsection). 

Consider now a continuous sheet DM with an unconventional influence function that is much 
narrower than the actuator spacings (green). Specifically, suppose it is a factor of M smaller than 
the blue influence function (in area for 2D). Its Fourier Transform intensity |-F{/}| 2 (Figure [8] 
middle right, green) shows its ability to correct in super-Nyquist. It will be a factor of M wider (in 
area for 2D) and a factor of M dimmer than the case of a conventional DM (blue), so as compared 
to a conventional DM, it can control super-Nyquist regions, at the expense of lowering the contrast 
of errors it can control in the sub-Nyquist region. Roughly speaking, such a DM is capable of 
correcting any one of M super-Nyquist regions, as well as the sub-Nyquist region, as long as the 
contrast of the errors is not higher than 1/MN, averaged across any one region. Note that we can 
create such an ’’unconventional” DM using existing conventional DMs. For example a conventional 
64x64 DM with only every 8th row and column connected is effectively an 8x8 DM with undersized 
influence functions. In this case, N= 8x8 and M= 8x8, the DM can still control speckles out to 32 
A/D, but only one 8x8 A/D-sized region at a time. If the number of electrical lines is a cost or risk 
driver, and a mission emphasis is on planet characterization rather than search (i.e. planet location 
is known), such a DM may be preferred to the fully connected 64x64 DM. 

Now consider the case of a conventional influence function (blue), but with a dip in the middle 
from print-through manufacturing (blue dotted curve in Figure [8]). This influence function was 
created simply by taking the difference between the blue and green influence functions. As a result, 
the corresponding (^\AEf\ 2 '/ has essentially the same controllability of the sub-Nyquist region as 
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the solid blue curve, but also can control super-Nyquist regions as well as the green case. The 
relationship between the width of the narrow dip in the influence function, the contrast and the 
number of super-Nyquist regions that can be controlled is the same as for the green case. 

Finally, a segmented DM (red color) has both a full control of the sub-Nyquist region (to full 
contrast levels similar to the conventional continuous-sheet DM case in blue) and can correct super- 
Nyquist regions, as long as the error does not exceed the contrast shown by the red curve when 
spatially averaged across any super-Nyquist region. This ability comes from the side lobes of the 
red curve. Correction will be best at odd multiples of half-Nyquist frequency (peaks of sidelobes) 
and there will be a blind spot at multiples of Nyquist frequency (zero-crossings of the red curve). 
In 2 dimensions, these sidelobes will be strong only along directions normal to the atuator edges 
(along ’’diffraction spikes” of the PSF). 


4.3. SNWC with a grating or beamsplitter 

DM diffraction into super-Nyquist regions can be caused not only by a particular shape of the 
influence function, but also directly by a mild external grating. The case of a grating with periodicity 
matching the DM actuator periodicity has already been treated in the previous subsection (dotted 
blue case in Figure [8]), and here we consider the case where the grating has a finer periodicity, 
which can be designed for a particular desired super-Nyquist distance and contrast level. 

In the previous subsection, the action of the diffraction-causing agent (influence function) was 
a convolution in the pupil plane and a multiplication in the image plane. In this subsection, the 
action of the diffraction-causing agent (grating) is the opposite: a multiplication in the pupil plane 
and a convolution in the image plane. 

Assume the same scenario as in the previous subsection, except with a mild grating represented 
by a periodic function g(x) multiplying the DM field ©• Define the new DM field as: 

Edm, 9 (x ) = E DM (x)g{x) (10) 

= A(x)(l + AE DM (x))g(x) 

= A(x)g(x) + AE DM (x)g(x) 

where E^m is the field from the previous subsection (i.e. without the presence of the grating). 
In the image plane, the first term will lead to the on-axis star PSF, together with fixed (DM- 
independent) attenuated copies of the star PSF diffracted by the grating. Each of these PSF copies 
can be thought of as a new star around which conventional sub-Nyquist wavefront control can be 
applied, thus enabling super-Nyquist wavefront control with respect to the original star. However, 
each PSF copy will create a small ’’blind spot” that the DM will not be able to remove, similar 
to the blind spot created by the original star in the center of the image. (If the periodicity of the 
grating is the same as the DM pitch, as would be the case if the grating is the DM print-through 
pattern, then there will be a small blind spot in the exact middle of every super-Nyquist region.) 
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For purposes of this subsection, we will assume that any light these PSF copies create outside the 
blind spots are part of the star leak error to be suppressed. Because they are independent of the 
DM, we do not book-keep them as part of the DM correction field and thus ignore the first term, 
focusing only on terms that are DM-dependent. The perturbation of the DM electric field by the 
DM is then: 


A E DM , g {x) = A E DM (x)g{x) (11) 

Fourier transforming to the image plane leads to the perturbation by the DM of the image plane 
electric field: 


a E f , g (0 = A D/(£)G(£) (12) 

where £ is the image plane position, A Ef is the image field perturbation from the previous sub¬ 
section (i.e. without the grating), and G = F{g}- Because g is a periodic function, its Fourier 
Transform has the form: G = )T) g n 5 (£ — nA£) where g n are the Fourier Series coefficients of g{x) 
and A£ is the spacing between the diffraction orders in the focal plane (equal to the spacing between 
super-Nyquist regions if the grating periodicity is the same as the DM actuator periodicity). 

Following the same statistical characterization method as in the previous subsection, we treat 
the perturbations of the image plane electric field by the DM as a random variable and compute 
its power spectral density, which after some algebra becomes: 

<|AD / J 2 ) = IV|.F{/}| 2 *|G| 2 (13) 

In other words, the sub-Nyquist control region is copied to the diffraction orders of the grating g 
and attenuated by a factor equal to the contrast of the diffraction order. Two examples are shown 
in Figure [9]. The blue case corresponds to the case where G only has one off-axis term, which is 
simply a mild beamsplitter bleeding off 1% of the total beam into 20 A/D. This is the most efficient 
way to control a specific fixed super-Nyquist region. The green case in Figure [9] corresponds to a 
periodic amplitude grating such as a print-through on the DM. It acts in much the same way as a 
beamsplitter, except it creates many super-Nyquist control regions instead of just one. 

There are a few key principles relating the grating design characteristics and the super-Nyquist 
control contrast and region location: (a) when a grating creates a diffraction order (PSF copy) of 
contrast C, super-Nyquist control is enabled around that diffraction order and up to IV/2A/D away, 
in exactly the same fashion as sub-Nyquist control around the original on-axis star; (b) the total 
energy of the error to be corrected cannot exceed the energy in that diffraction order. For example, 
if the diffraction order is 10 -3 contrast and we have a 32x32 DM correcting in a 32x16 A/D half¬ 
region around that diffraction order, then the average contrast of the error corrected cannot exceed 
1CD 3 /(32x16) = 2 x 1CP 6 . (On the other hand, if the region of interest is only 3x3 A/D, then 
speckles of up to 10 -4 contrast can be corrected. 
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Fig. 9.— Grating case. The left column shows the pupil plane and the right column the image 
plane. The top row shows the DM shape and its Fourier Transform (Fourier series of a n ). The 
second row shows the grid pattern and the bottom row shows the multiplication or convolution of 
the top two rows. 


4-3.1. SNWC implementation with a grating 

Using dots on the pupil has been proposed previously to calibrate dynamic distortions on wide- 
held optical systems enabling high-precision astrometric measurements (Guyon et al. 2012). For 
this technique the dots can be arranged in the pupil using a hexagonal geometry allowing higher 
azimuthal sampling. The diffractive pupil spacing can be adjusted to create PSF replicas to run 
the SNWC, and also to obtain high-precision astrometry on wide-held images. A description of the 


optimal hexagonal geometry for combining these two techniques has been published by (Bendek 


et al. 2013). 
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4.4. Multi-Star Wavefront Control 


For sake of completeness, we will describe here how to simultaneously suppress the speckle 
fields of both stars (MSWC). Consider the case of two stars, A and B, with (for now) a sub- 
Nyquist separation. The main challenge is that light from the two stars is mutually incoherent, 
and therefore light from each star can only be used to destructively interfere its own speckle field 
but not the field from the other star. In order to suppress both stars, it is necessary to be able 
to independently modulate the speckle field of each star without affecting the speckle field of the 


other (in some region of interest). Figure 10 shows the special case of two stars separated by 16 
A/D. The left panel shows the (sub-Nyquist) control region of a 32x32 DM with respect to star 
A, and separates the region into 4 vertical sections, each controlled by a different and independent 
set of modes on the DM (the outer sections are controlled by modes on the DM corresponding to 
spatial frequencies of 8-16 cycles per aperture (cpa), and the inner regions are controlled by 0-8 


cpa). The middle panel of Figure 10 shows the same thing, but with respect to star B. Finally, the 
right panel superimposes these control regions of the two stars. The two regions between the two 
stars are labeled as I and II . In region I, the 0-8 cpa modes modulate the speckle field of the B star 
but not the A star, and the 8-16 cpa modes modulate the speckle field of the B star but not the A 
star. In other words, in region I, the speckle field of the A star can be suppressed by using the 8-16 
cpa modes without affecting the B star and the speckle field of the B star can be independently 
suppressed by using the 0-8 cpa modes without affecting the A star. In effect, we have reduced 
the MSWC problem to two separate conventional WC problems (each using different DM modes), 
which we know how to solve. Solving these two conventional WC problems simultaneously results 
in simultaneous suppression of the speckle fields of both stars in region I. The same can be done 
in region II (but not simultaneously with region I). The final result is that a double star dark zone 
can be created at the expense of reducing the size of the control region by a factor of two (once 
again, one cannot cheat the number of degrees of freedom available on the DM), and there are two 
such regions. These two smaller regions can be suppressed separately and then stitched together 
to create the full field of view between the two stars. 

This idea can be generalized to the case of arbitrary (sub-Nyquist) star separation. The result 
is that one can always partition the intersection of the sub-Nyquist regions of the two stars into 
two sections, each with half the area of the original control region, where speckle fields of both stars 


can be simultaneously suppressed. (These regions will have shapes different from Figure 10 and 


may consist of disconnected parts.) A generalization to N stars, implies N independent correction 
regions, each with 1/N of the original control region area. These regions are found to be so-called 
Voronoi partitions of a periodically extended star field. It should be noted that in practice it is 
impossible to completely decouple DM modes. Any DM mode will always affect all stars everywhere 
to some level. However, on the regions we constructed, one of the stars speckle fields is affected 
much more than the other, so the above algorithm works in closed loop. 
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Fig. 10.— Use of different DM modes to independently modulate light for two different stars. The 
first scheme on the left shows the different controllable regions for the star A (black dot) and the 
middle the controllable regions for the star B (gray dot). The figure on the right represents the 
superposition of the two and shows that there are two regions names I and II between star A and 
star B for which we can control A and B at the same time with the same DM. For region I we 
would use the modes from 0 to 8 A/D for B and 8 to 16 X/D for A and vice et versa. 


4.5. Super-Nyquist Multi-Star Wavefront Control 


In this section, we combined SNWC and MSWC and treat the case of two (or more) stars 


having a super-Nyquist separation. Figure 11 shows a diagram of an on-axis star A (suppressed 
by a coronagraph if present), and the companion off-axis star B. Both stars have sub-Nyquist 
control regions where conventional WC can suppress the speckles field of one star but not the 
other. Suppose that the mild DM grating diffracts the (attenuated) replica of star B inside the 
sub-Nyquist region of star A, just as in the case of SNWC. As discussed in section |4.1[ we can 


treat this diffracted replica as an actual star, effectively resulting in two stars with a sub-Nyquist 
separation. This reduces the problem to that of MSWC (which can be rigorously formalized at least 


in monochromatic light), which we demonstrated how to solve in section 4.4 The main potential 
difficulty arises for the case of broadband light, where the diffracted replica of the off-axis star looks 
spectrally elongated. This will certainly reduce the size of the correction region, but if SNWC works 
in broadband light then SNMSWC also works in principle with a reduced dark zone size. 
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Fig. 11.— Diagram of Super-Nyquist Multi-Star Wavefront Control. The DM grating diffracts an 
attenuated replica of star B into a sub-Nyquist region of star A, just as in the case of pure Super- 
Nyquist Wavefront Control. This brings the problem back to the simple Multi-Star Wavefront 
Control, for which we treat the diffracted replica as an actual star. In this image a coronagraph 
blocks the light originating from star A. A side effect that can be seen on the diagram is the replica 
of A in the controllable region of B. This allows us to then search for planets around A. 


5. Simulation of SNWC 

To demonstrate the SNWC, we chose to simulate the observation of binary of equal brightness 
separated by 100 X/D. This could be the components of a Cen with a 1.5m telescope (such as Exo- 
S Q) in monochromatic as well as polychromatic light (10% bandwidth with a central wavelength 
of 770nm). The expected separation of a Cen components A and B in the year 2025 is about 10 
arcsec, which corresponds to about 100 X/d away at a wavelength of 770 nm. This is also applicable 
to imaging large disks. Since this paper demonstrates SNWC and not MSWC, we are assuming 
that the light from the target on-axis star has been completely removed (e.g. by a starshade) and 
only the off-axis companion remains. The goal is to demonstrate in simulation that it is indeed 
possible to create a dark zone beyond the Nyquist frequency of a DM in this configuration. 


2 http://exep. jpl.nasa.gov/stdt/exos/ 
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5.1. Simulation description 


In this paper, we have focused on demonstrating wavefront correction, assuming we have 
perfect knowledge of the phase and amplitude of each star’s speckle field in the region of interest. 

The problem of correction is conceptually separate from the problem of estimation and is in 
some sense more fundamental, because the ability to correct implies estimation is possible, but the 
converse does not hold. 


Conventional WC estimation (e.g. in EFC or stroke minimization, (Give’on et ahj|2007; Pueyo 


et al. 2010)) uses the DM to create known probe fields to modulate the region of interest and 
then analyzes this modulation to reconstruct the phase and amplitude of the pre-existing speckle 
field. Similarly to the correction problem, the estimation problem in all cases (SNWC, MSWC, 
SNMSWC) can also be reduced to conventional wavefront estimation as long as the DM is able to 
create known probe fields in the region of interest of sufficient energy, which will be the case as 
long as we demonstrate that correction is possible. 


In section 4.1, we described that SNWC requires the DM to diffract light beyond the Nyquist 
region with, for example, the periodic print-through pattern left over from manufacturing. To create 
such a pattern, we used a grid of dots in the image plane, created from a mask of uniform intensity 
equal to 1 with a grid of 0. The period of the grid is set such that there will be a diffracted (and 
attenuated) PSF copy next to the zone we would like to correct. However, even with a fixed number 
of actuators and therefore a fixed frequency created, (multiple of the number of actuators (here 
32)), we can cover all separation scenarios. Another parameter that can be adjusted in order to 
control the performance is the width of the dark zone region. Indeed the bigger the region, the more 
DM stroke is needed to achieve deeper contrast. We used a 4x8 \/d region for the demonstration, 
which we found being a good compromise between the performance and the discovery region. The 
simulations were done both in monochromatic light (770 nm) and in polychromatic light with a 
10% band. Finally, we also studied the case of a non-aberrated and aberrated wavefront. The 
aberrations were introduced in the pupil plane as a power law with a coefficient equal to -2. 


5.2. Results 


Monochromatic Light In monochromatic light, the introduced diffraction grid creates dots at 
100 A/d with an intensity of 1.34xl0~ 3 relative to the central star. Figure 


12 


shows the results 

of a (nonlinear) correction solution without aberrations (left) and with 25nm rms of aberrations. 
The median contrast obtained without any aberrations is 2.5xlO~ 10 and with aberration we reach 
a contrast of 1.7xl0~ 9 . This corresponds to a factor 100 improvement from the no grid simulation 
for both the aberrated and non-aberrated case, demonstrating the ability of SNWC to create dark 
zones outside the nominal Nyquist limit of the DM. We expect that deeper contrasts are possible 
with better nonlinear solutions, but such optimizations are outside the scope of this paper. 
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Fig. 12.— Super-Nyquist simulation results in the monochromatic case, for a region of interest of 
4x8 X/D, Left: Onm rms aberrations. Right: 25nm rms aberrations. The median contrast obtained 
is 2.5xlO -10 without any aberrations and 1.7xl0~ 9 with aberrations. 


Polychromatic Light We now consider a more realistic scenario and study the effect of poly¬ 
chromatic light. We use the methods in (Give’on et al. 2007) for polychromatic control within a 
10% bandwidth. In order to get an accurate image, we chose to sample the 10% bandwidth at 3 
wavelengths. This allows a good compromise between computational speed and spectral resolution. 
For better contrasts or larger bandwidth, one will need to increase the sampling of the bandwidth. 
Figure 13 shows the results without and with lOnm rms aberrations. The median contrast ob¬ 
tained without any aberrations is 4.9 x lO -9 and very similarly for the case with aberrations, we 
reach 5.3xlCP 9 . This gives the hint that we are limited by chromaticity. Keep in mind however 
that these results are raw contrasts before any type of post-processing techniques that would allow 
us to gain a factor 10 or even 100 on the contrast, which is encouraging for planet detection and 
characterization. 


Table [l] shows the summary of the different simulation results in the case of a equal brightness 
system. 

This table shows raw contrast, before any post-processing algorithms. Therefore, it is very 
encouraging to see that it would be to detect Earth-like planets around both components of a 
Centauri. 
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Fig. 13.— Super-Nyquist simulation results in the polychromatic case, 10% bandwidth around 
550nm. The dark zone region is around 4 and 8 A /D. The figures on the left show the case of 
no aberrations and the figures on the right with lOnm rms of aberration. The median contrast 
obtained is 4.9x10-9 without aberrations and 5.3x10-9 with aberrations. This shows we are limited 
by polychromaticity. 


6. Conclusion 

In this paper, we presented a method that will enable high-contrast imaging of multiple stars 
systems. This new concept was broken in in three different steps, namely; (a) Super-Nyquist 
Wavefront Control (SNWC); (b) Multi-Star Wavefront Control (MSWC); (c) Super-Nyquist Multi- 
Star Wavefront Control (MSSNWC). We showed with simulations that it is possible to create a 
dark zone past the Nyquist frequency of a deformable mirror using a diffractive grid in the pupil. 
This grid can be the DM itself or an additional mask. Performance was shown in the case of a 
Cen, an interesting target since it is the closest potential earth-like planet hosts. This particular 
scenario can be applicable to a perfect star shade coronagraph for which we totally block the light 
coming from the parent star but still need to remove the light coming from the companion in the 
dark zone of interest. The next step is to show that the DM can create a dark zone when adding 
both components in the lab. This will be done either using one or two DMs. The main scientific 
impact of this work is to enable direct imaging of planetary systems and disks around multiple star 
systems as well as in regions far from the star. This can be done at little additional hardware cost 
or changes to existing mission concepts, such as AFTA, Exo-C, Exo-S, and EXCEDE will greatly 
multiply the science yield of these missions. 


The material is based upon work supported by the National Aeronautics and Space Adminis¬ 
tration under Prime Contract Number NAS2-03144 awarded to the University of California, Santa 
Cruz, University Affiliated Research Center. This work was supported in part by the National 
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Table 1: Summary of performance achieved with the simulations before and after wavefront control 
(WFC) and for monochromatic and polychromatic cases. 

Aberration Contrast 


No Grid, Monochromatic, before WFC 

Onm 

3.5 

10- 

-8 

No Grid, Polychromatic(10%), before WFC 

Onm 

3.6 

io- 

-8 

Grid, Monochromatic, after WFC 

Onm 

2.6 

10“ 

10 

Grid, Polychromatic(10%), after WFC 

Onm 

4.9 

io- 

-9 

No Grid, Monochromatic, before WFC 

25 nrn 

5.4 

io- 

-7 

No Grid, Polychromatic (10%), before WFC 
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8.7 
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-7 

Grid, Monochromatic, after WFC 
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4.8 

io- 

-9 
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Aeronautics and Space Administration’s Ames Research Center under a center investment fund. It 
was carried out at the NASA Ames Research Center. Any opinions, findings, and conclusions or 
recommendations expressed in this article are those of the authors and do not necessarily reflect 
the views of the National Aeronautics and Space Administration. 
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